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ABSTRACT 

We investigate which practical constraints are imposed by foregrounds to the de- 
tection of the B-mode polarization generated by gravitational waves in the case of 
experiments of the type currently being planned. Because the B-mode signal is prob- 
ably dominated by foregrounds at all frequencies, the detection of the cosmological 
component depends drastically on our ability for removing foregrounds. We provide 
an analytical expression to estimate the level of the residual polarization for Galactic 
foregrounds, according to the method employed for their subtraction. We interpret 
this result in terms of the lower limit of the tensor-to-scalar ratio r that allows to 
disentangle the cosmological B-mode polarization from the foregrounds contribution. 
Polarized emission from extragalactic radio sources and gravitational lensing is also 
taken into account. As a first approach, we consider the ideal limit of an instrumental 
noise-free experiment: for a full-sky coverage and a degree resolution, we obtain a 
limit of r 10~^. This value can be improved by high-resolution experiments and, 
in principle, no clear fundamental limit on the detectability of gravitational waves 
polarization is found. Our analysis is also applied to planned or hypothetical future 
polarization experiments, taking into account expected noise levels. 

Key words: cosmic microwave background — polarization — cosmological param- 
eters — early Universe 



1 INTRODUCTION 

Standard cosmology is based on the inflationary paradigm, 
which provides a good description of the observed Universe. 
Besides a flat Universe, inflation predicts a Gaussian, adia- 
batic and nearly scale-invariant spectrum of primordial per- 
turbations. However, inflation is lacking a clear physical sup- 
port and alternative theories leading to the prese nt Universe 
have been investigated (e.g.. iKhourv et alj|200lf) . 

A further prediction of inflationary models is the pres- 
ence of tensor perturbations to the spatial metric (gravita- 
tional waves). The detection of these perturbations would 
be an extraordinary way to distinguish inflation from other 
competing scenarios and, at the same ti me, to constrain 
the p roperties of the inflationary potential jPodelson et alJ 

The Cosmic Microwave Background (CMB) 
anisotropies and polarization probably offers the best 
means (if not the only) of detecting the tensor metric 
perturbations. In fact, in most inflationary models, the 
amplitude of the power spectrum for "tensor" temperature 
anisotropies is directly related to the energy scale at 



which the inflation has occurred. Considering a single-field 
infiation (to lowest order in the slow-roll parameters), the 
relation between the infi ationary potential, U(( /)), and the 
tensor quadrupole is iTurner fc Whit3ll99d) 

■^^1.65Ql, (1) 

■mpi 

where U* is the inflationary potential evaluated when the 
present Hubble scale crossed outside the horizon during in- 
flation. Then, the "energy scale" of the inflation Ei — Vt^^ 
can be expressed in terms of the scalar quadrupole Q| and 
the ratio of tensor to scalar quadrupoles r — Qt/Q%^ as 

~ Q^/'r'/^mp, ~ 3 X lo" r'^^ GeV , (2) 

where we have used the COBE estimate Qs — 18jj,K 



^ The tensor-to-scalar ratio is often defined also in terms 
of the primordial amplitude of tensor and scalar fluctuations 
iLeach et al.ll2002ft . For the "concordance" model, it corresponds 
to 1.78Qf,/Q|. 
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iHinshaw et aljfl99(J) . However, tensor perturbations con- 
tribute to the temperature power spectrum essentially at 
low I, with a roughly scale-invariant spectrum. Unless the 
relative amplitude of tensor perturbations respect to scalar 
ones is large enough (r > 0.1), measurements of the tem- 
perature power spectrum can only provide an upper limit 
on due to constra ints imposed by the cosmic variance 
f Knox & Turner 1994) [for example, WMAP finds r < 0.71 
tSpergel c t al. 200^)]. 

The detectability of the contribution of tensor per- 
turbations is much more promising exploiting CMB po- 
larization measurements. The polarization field can be di- 
vided into a curl-free component of even parity called 
"E-mod e", and a curl component of odd parity, "B- 
mode " jKamionkowski et alJ Il997l : IZaldarriaga fc SeliakI 
^93). While the former is generated by both scalar and 
tensor perturbations, B-mode polarization arises only from 
tensor perturbations and carries thus a direct imprint of the 
inflationary epoch. Nevertheless, the B-mode intensity is ex- 
pected to be extremely weak, and, even for optimistic values 
of r, the rms signal is only a fraction of /xK, less than 1% of 
the level of temperature anisotropics at degree scales. Future 
experiments will need to reach extraordinary sensitivities to 
be able to detect such a low signal. 

Gravitational waves produced during inflation are not 
the only possible source of primordial B-mode polariza- 
tion. For example, the existence of tangled primordial 
magnetic flelds or cosmic strings would generate a vector 
component of metric perturbations, providing a contribu- 
tion to the B-mode polarization on small angular scales 
(^Seshadri fc Subr amanianI l200ll : ISubramanian et alJ 120031: 
FPoeosian et aLi2003h . In anv case, the discussion of these 
sources of polarization is outside the scope of this work. 

Even if the instrumental sensitivity appears to be the 
biggest problem for detecting the tensor-induced polar- 
ization, important constraints come also from the pres- 
ence of other sources of the B-mode polarization, as 
foregrounds and effects mixing CMB E- and B-modes. 
The most relevant mode-mixing effect is due to the 
gravitational lensing produced by large scale structures, 
which conve rts a fraction of the CMB E -mode component 
to B-mode JZaldarriaga fc Seliak^ll998^. Different authors 
JHu fc Okamotoll200l iKesden et alJl2003l ISeliak fc Hiratal 



l2003l) presented methods to reconstruct the gravitational 
lensing potential using information from the CMB polar- 
ization itself and to remove the lensing contamination. The 
problem of a correct separation between E- and B-mode 
in the presence of a partial sky coverage, pixelization and 



ature (seelLewis et al.ll2002 


. lBunnll2002l iBunn et alJl2003l 


iHu et alJl2003l.lLewidl2003l. 


Challinor fc Chonll2004). and a 



suitable formalism to minimize their influence has been pro- 
vided. 

Another critical source of confusion for the detection 
of the CMB polarization is expected from foreground con- 
tamination. Apart from the free-free emission that is not 
polarized, foregrounds have been shown to have, in general, 
an high degree of polarization equally shared between E- 
and B-mode: synchrotron emission is on average 20%-30% 
polarized while dust and extragalactic sources are polarized 
from few per cents to 10 per cent or more. Foregrounds are 
expected to dominate the sky B-mode polarization at all 



frequencies and angular scales, even for high values of r. 
On the other hand, they are also a serious problem for the 
CMB E-mode signal which is at most a 10% of temperature 
anisotropics, but significantly less on large angular scales. 

A secondary B-mode polarization is also generated by 
scattering of the CMB photons from ionized gas in galaxy 
clusters or in a patchy reionizat ion sc enario (see lHull200d . 
iLiu et alJ 1200 ll iBaumann et alJ [2003', San tos et al.l l200:j^. 
However, these effects produce a polarization that is several 
orders of magnitude below the dominant contributions. 

The main question that this paper address is how Galac- 
tic and extragalactic foregrounds (including the gravita- 
tional lensing contamination) limit the detectability of the 
CMB B-mode polarization. In terms of cosmological con- 
straints, we want to estimate the lower limit of r (hence of 
the energy scale of infiation) under which the signal from 
tensor perturbations cannot be distinguished by polariza- 
tion measurements. The case of an ideal experiment (i.e., 
noise-free experiment) is first considered. Then, the same 
analysis is also applied to planned or hypothetical experi- 
ments, giving hints on the best strategy for future missions. 



2 INTENSITY LEVEL OF B-MODE 
POLARIZATION FOR CMB AND 
FOREGROUNDS 

In this section we discuss the expected level of the polarized 
signal for the CMB and foreground emission as a function 
of the observing frequency u and the angular scale. The re- 
sults are presented in Figures and |21 The left plot of Fig- 
ure0shows the dominant foregrounds in the l-u space. The 
not-shaded area is where the CMB E-mode is higher than 
the foregrounds E-mode contribution. The diffuse Galac- 
tic emission (synchrotron and dust emission) is essentially 
the strongest foreground in polarization. Only at high I and 
frequencies lower than 100 GHz extragalactic radio sources 
start to be relevant. We notice how the "cosmological win- 
dow" (the l-u region where CMB exceeds foregrounds emis- 
sion) is strongly reduced as compared to temperature fiuc- 
tuations, expecially at large angular scales. In the right plot, 
we consider the polarization level of Galactic emissions but 
reduced by an order of magnitude (this is, at least, what we 
expect to reduce in a multyfrequency experiment). In this 
case, extragalactic radio sources become an important con- 
tribution already for I ^ 100, as well as the polarization in- 
duced by gravitational lensing at 100 ^ v ^ 200 GHz. Now 
the not-shaded area represents the £-v region dominated 
by the CMB B-mode polarization for a cosmological model 
with r = 0.1. As expected, it is confined only to ^ ^ 100. 

In Figure |5| we plot the rms value for the X = E, B 
modes polarization^ as expected to be measured by an ex- 
periment with a resolution of 1°: 



+ i)Cxiw! 



1/2 



(3) 



where Wi is the beam function of a hypothetical instrument 
and Cxi is the angular power spectrum (APS) of CMB or 



^ Hereafter, we are supposing that, on average, the amplitude of 
E- and B-mode polarization is the same for foregrounds emission. 
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foregrounds. The rms values are given in antenna tempera- 
ture. The 1° resolution is chosen to estimate the rms polar- 
ization because degree angular scales are the most interest- 
ing for detecting the CMB B-mode polarization. Let us see 
now the properties for each polarized component in detail. 



2.1 Cosmic microwave background 

The CMB B-mode power spectrum (Cbi) is characterized 
by a peak at Ip ~ 90, i.e. at the angular scale corresponding 
to the horizon at recombination. The amplitude of the spec- 
trum is directly proportional to the cosmological parameter 
r and is related to the energy scale of inflation: at the peak 
the amplitude is 



ip(tp + 1) 
27r 



Cbi 



1/2 



To 



0.3 /^Vk 0.03^ 



1016 GeV 



(4) 



At I > Ip the spectrum rapidly decreases because grav- 
itational waves oscillate and decay once inside the hori- 
zon. At very large scales very low B-mode polarization is 
also expected because, at the moment of recombination, the 
anisotropy quadrupole moment has not been significantly 
produced yet on scales larger than the horizon. However, 
the reionization of the Universe has the effect to partially 
polarize the CMB radiation and shift to lower redshift a 
fraction of the last scattering, producing a further peak at 
low I in the spectrum. The amplitude and the position of 
such a peak depends on the opti cal depth of the Un iverse, 
r [WMAP finds r = 0.17 ± 0.07 iSpergel et al.ll2003M . 

In Figure 2b we report the expected values of Brms for 
r = 0.1, 10~^ and 10~^ (solid red lines). As shown in Eq. 
(0J|, it scales like r^^'^. If only very large scales are consid- 
ered, -Brms is strongly dependent on r and exceeds 0.1 
for r and r > 0.1. On the contrary, at degree resolution, r is 
less relevant but not negligible (-Brms increases by a factor 2 
changing r from to 0.17). In the plot we consider r = 0.1 
(a conservative value compared to the WMAP best-fit one). 
The B-mode spectrum is less sensitive to other cosmologi- 
cal parameters. We fix those in agreement with the "concor- 
dance model" (Slo = 1, = 0.7, Q.b = 0.05 and h = 0.7). 

The detection of the CMB E-mode polarization is much 
easier than that of B-mode and it is even better with sub- 
degree resolution experiments (see Figure [P. At these angu- 
lar scales, DASI llKovac et alJl2002^ achieved the first direct 
measure of E-mode polarization, while more recent obser - 
vations hav e been provided by DASI jLeitch e t al.ll2004). 
CAP MAP jBarkats et alJl2005^ and CBI llReadhead et alJ 
|20^).At degree angular scales the value of -Erms is less than 
1/iK with a polarization degree of only few per cent or less. 
Nevertheless, the CMB E-mode is still the dominant con- 
tribution in a wide range of frequencies, approximately be- 
tween 60 < < 150 GHz. 

Finally, we plot in Figure 2a the level of the B-mode 
polarization induced by the gravitational weak lensing of 
large scale structures (see magenta line). This estimate is 
obtained using the CMBFAST package^. 



http: / /www. cmbfast.org 



2.2 Synchrotron emission 

Among all foregrounds, synchrotron emission is the most 
polarized (up to ~ 75%). It is expected to be the dominant 
component at low frequencies ^ 70 GHz). Total-intensity 
data are available up to ~ 20 GHz, while polarization obser- 
vations of the synchrotron are obtained only at few GHz 
or less. Estimates on the amplitude of the polarized syn- 
chrotron emission at cosmological frequencies are especially 
complicated for two reasons: the steepening of the spectral 
index that may occur ai v ^ 20 GHz, and the frequency 
dependence of the degree of polarization. 

The synchrotron spectral index, /3s, is computed 
on nearly all sk y using th e surveys at 408 MHz 
jHasIam et alJll982^. at 1.4G Hz fReich fc Reicblll986[l and 
at 2.3 GHz iJonas et ai]|l998h bv lCiardino et al.N2002f) and 
iPlatania et al.l (l200^i ~These analysis show Ps to vary across 
the sky between —2.5 and —3.0, with an average value < 
/3s >~ —2.7 and a dispersion of 0.1-0.2 at sub-degree angu- 
lar resolution. A relevant steepening (< /3s >~ —3 out of the 
Galactic Plane) is found from the first- year dat a of WMAP 
at 23 GHz ( Bennett ct al. 2003). Bcrnar di et al.l (|m04) stud- 
ied the distribution of the spectral index computed be- 
tween 1.4 and 23 GHz: the frequency spectrum is observed 
rather fiat along the Galactic Plane, whereas steeper at high 
Galactic latitudes with a typical value /3s = —3.1. A fur- 
ther stccpcniii g is al so suggested in the 23-41 GHz range by 
[Bennett et alJll2003l) . although it has to be confirmed. In the 
following, we shall assume /3s — —3.1 as the mean spectral 
index for the synchrotron emission at all considered frequen- 
cies. 

Estimates of the rms signal oi X = E, B modes for the 
synchrotron emission are reported in Figure 2a (green lines) . 
A set of estimates obtained in different ways is considered, 
which permits to get a range of plausible levels. First we 
scale down unpolarized synchotron as follow: the full-sky 
synchrotron map at 23 GHz, obtained from WMAP data 
after co mponent separation with the Maximum Entropy 
Method teenn ett et ai]l2003l) and smoothed by a Gaussian 
beam with FWHM= 1°, is used to compute the rms am- 
plitude of temperature fluctuations at high latitudes (|b| > 
20°). We find AT^ms = 82.6. From the iBrouw fc Spoelstral 
(|l973) observations (the only up-to-now available data of 
the synchrotron polarization on a large fraction of the sky) , 
Spoclstra ( 1984) found that on degree scales the fractional 
polarization at 1.411 GHz is typically 10-20 per cent with 
the highest value being 35%. Taking into account that Fara- 
day rotation can reduce the polarization degree at this fre- 
quency, we assume that synchrotron emission is polarized 
between a 10 and a 30 per cent as lower or upper limit. We 
get estimates as plotted with dotted lines in Figure 2a. 

Secondly, estimates of Xr ms can be direc tl y ach ieved 
also from polariza t ion d ata. IBruscoli et al.l (|20o3) use 
iBrouw fc Spoelstral lll97d) survey to compute the E- and 
B-mode APS in three patches with high signal-to-noise ra- 
tio centered at the latitudes 6 = 5°, 44.5°, 72.5°. They find 
that the APS can be well approximated by a power law with 
a quite fiat slope, ^ —2. Taking as the "average" spectrum 
at 1.4 GHz Cxe = Omr^-^{K^), we find Xrin3(1.4GHz) = 
O.llK at 1° resolution (dashed line in Figure 2a). Sim- 
ilar results are also obtained from the ana lysis of high- 
resolution surveys at low Galactic latitudes (iDuncan et alJ 
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Figure 1. The dominant foreground in polarization as function of the frequency and the angular scale (synchrotron: light blue; dust: 
yellow; extragalactic radio sources: brown; gravitational lensing: green). Galactic foregrounds are modelled by Eq. I14i - I15i (see also 
Figure 2b). For extragalactic radio sources we use a flux limit of 1 Jy. Left plot: the not-shaded area is where the CMB E-mode is the 
dominant polarized signal. Right plot: as in the left plot but after reducing by a factor 10 the synchrotron and dust emission. Now, the 
not-shaded area indicates the region where the CMB B-mode (considering r = 0.1) is higher than the foregrounds contribution. 
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Figure 2. The rms value of foreground and CMB X=E, B modes vs. the frequency, as expected by an experiment with FWHM= 1°. 
—a- Estimates of foregrounds from different data. Synchrotron (green lines): 10%-30% of the WMAP ATrms (dotted lines); estimates 
from Brouwn & Spoelstra dat a (dashed line), fro m high-resolution low-latitude polarization surveys (solid line) and from observations 
of high-latitude polarization jAbidin et al ]|2003f: long-dashcd line). The spectral index is /3s = —3.1. Dust (dark blue lines): 5% of 
the WMAP ATrms (dashed line) and the result from the Prunct ct al. ( 1998) model (solid line). For the frequency spectr um we consider 
a one — component dust model with T = 18 K and emissivity a = 1.7. The dotted line is obtained by the "100/^m map" I' Schlegel et alj 
Il998ll extrapolated to m icrowave frequencies using the model 8 of| Finkbeiner et a l." 1 1999). Radio sources (light blue lines): estimates 
from lTucci et aTl l200'd) . taking the flux limit 1 and 0.2 Jy. Lensing— induced polarization (magenta line). 

— b— The black solid line is the total Xrms c ontribution of for e groun ds, taking for the synchrotron the 20% of the WMAP estimate 
and /3s = —3.1; for the dust the model 8 of iFinkbeiner et al.l ll999l) . CMB (red lines): Erms (dotted line) and Brms (with r = 
0.1, 0.01, 0.001; solid lines) for the "concordance" model with r = 0.1. 
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119971. iDuncan et ai]ll999l lUvaniker et"aIlll999^: extr a polat - 
ing the polarization AP S provided by iTucci et al.l (|200(J) 
and lBruscoli et alJ lj200 2^ to low ^s, we find Xrms ~ 0.02 K 
at 2.4 GHz (solid line in Figure 2a). However, these esti- 
mates probably have to be considered as upper limits be- 
cause they come from observations of the Galactic Plane or 
of highly-polarized regions. Their very good agreement with 
the 30%-WMAP predictions shown by Figure 2a confirms 
that a polarization degree up to 30% can be observed in 
highly-polarized area but as it may be not the typical per- 
centage on the sky. The upper lines in Figure 2a may be 
somewhat pessimistic estimates. 

Finally, the recent observations of the high-Galactic lat- 
itude polarization w ith the EfTelsberg Telescope at 1.4 GHz 
jAbidin et alJl2003ft point out structures in polarization on 
scales of several degrees with a fractional polarization up to 
30-40% in the brightest regions and a rms Q and U signal 
of 50 mK. This value, extrapolated to WMAP frequencies, 
corresponds to the WMAP result with an average polar- 
ization degree of 10% (see long-dashed line). Nevertheless, 
as the polarization level at 1.4 GHz can be significantly re- 
duced by Faraday depolarization, a percentage of 20% for 
the polarized synchrotron emission may be more realistic at 
microwave frequencies. The lower lines in Figure 2a are thus 
probably somewhat optimistic. 

2.3 Dust emission 

Millimetric or submillimetric measurements of the Galactic 
dust polarization are usually concentrated in Galactic clouds 
and star -forming regions with arcminutes angular resolu- 
tion fe.g. lHildebrand et al.|[l99llGreaves et alJll999^ . They 
show a few per cent polarization and no clear frequency de- 
pendence. The first observations on large angular scales of 
the polarization of the G alactic dust emissi on are provided 
by Archeops at 353 GHz iBenoit et al.ll2003^ . The Archeops 
data show a significantly large-scale polarized emission in 
the Galactic Plane, with a polarization degree of 4-5%, with 
several clouds of few square degrees appearing to be polar- 
ized at more than 10%. 

Total-intensity observations at infrared wavelengths, 
where the thermal dust emission is very dominant, have 
been obtained by the I RAS and CQBE (DIRB E) satellites. 
Combining these data, ISchlegel et al.l il998h generated a 
full-sky map at 100/im with sub-degree resolution. A ten- 
tative extrapolation of the IQQtira map to microwave fre- 
quencies was done bv lFinkbeiner et alj (^9^ , assuming dif- 
ferent dust emissivity models. They use the FIRAS data 
in the 100-2100 GHz range to constrain dust properties. 
The best agreeement is met with a two-component model, 
consisting of a mixtur e of silicate and carbo n-dominated 
grams (see Model 8 of iFinkbeiner et al.llT99gl) . Extrapolat- 
ing the 100/.im map to 100 GHz with this model, we find 
ATrms = 7.4 fxK for a 1° resolution. The dust rms intensity 
at this resolution can also be estimat ed from WMAP data : 
using the Galactic emission models bv lBennett et all i2003l) , 
it is found ATrms = 13 /jK at 93.5 GHz, significantly higher 
than the Finkbc iner et alJ li 1999.1 value. However, because 
in the WMAP frequency range the dust emission is only a 
secondary contribution, this value is still controversial. The 
predictions on the rms polarization presented in Figure 2a 
assume a typical polarization degree for the interstellar dif- 



fuse dust emission of 5%, in agreement with Archeops results 
(dotted and dashed blue lines). 

Finally. fPrunet et alj (^9^ provide a model of the dust 
polarized emission using the three-dimensional HI maps of 
the Leiden-Dwingeloo survey at high Galactic latitudes. Dif- 
ferent hypothesis on dust grain properties and alignment are 
assumed. As main result, they estimate the dust polariza- 
tion APS at 100 GHz to be Cxt lO-^^"^ '' (/iK)^ corre- 
sponding to Xrms = 0.051 /^K at 1° resolution (solid line in 
Figure 2a). This value is significantly low compared to pre- 
vious ones and require a mean polarization degree less than 
1% for being consistent with total-intensity results. 

2.4 Extragalactic sources 

Extragalactic radio sources contribute significantly to tem- 
perature and polarization fiuctuations only at small angu- 
lar scales. At degree scales their polarized signal is a small 
fraction of the total foregrounds, even if it can exceed the 
cosmological B-mode for a low value of r. Predictions on 
the polarization properties of radio sources at microwave 
frequencies and on their contamination for CMB measure- 
ments are provided by Tucci ct al. ( 2004), using the evolu- 
tionary model of lToffolatti et alJ lll99lf) which predicts the 
number counts of sources at cm and mm wavelengths. In 
Figure 2a we report the rms polarization assuming that all 
sources with flux density higher than 1 or 0.2 Jy have been 
completely removed (light blue lines). 

At frequencies higher than 100 GHz, the contribution of 
dusty galaxies should also be taken into account. However, 
the sparsity of data on polarization at sub-mm wavelengths 
makes reliable predictions difficult. In any case, we expect 
that the polarization degree from the dust emission in exter- 
nal galaxies is not higher than the polarization degree ob- 
served in the Milky Way, i.e. only a few per cent. Supposing 
a polarization of 2% llHildebrandl fl996: Gre aves et ai]l2000l : 
iMatthews fc WilsonI [20021) . Figure ITI shows that the dusty 
galaxies contribution never overcomes the other foregrounds 
at all frequencies lower than 1000 GHz and at ^ < 3000. Po- 
larization from dusty galaxies will be neglected in this work. 

2.5 Other sources of polarization 

Here we give a brief description of some secondary sources 
of polarized emission, although they will not be taken into 
account in the analysis. 

The existence of an anomalous dust-correlated emis- 
sion, far brighter than the expected thermal d ust emission, 
has been well detected at ^10-30 GHz (e.g.. iKogut et alJ 
Il996llde Oliveira-Costa et alJll999l . iFinkbeinerlbOOSTl . One 
possible explanation is the electric dipole emission from 
rapidl y rotanting small dust g rains, known as "spinning 
dust" jPraine fc Lazarianlll998l) . This emission is expected 
to be polarized if grains are aligned with the magnetic field. 
In any case, as shown bv lLazarian fc Finkbeineiri2003t) . the 
polarization is marginal for > 35 GHz. Another possi- 
bility is the m agneto-dipole emission fr om strongly mag- 
netized grains dPraine fc LazariarJll99gl) . The polarization 
due to this mechanism depends on the dust composition and 
structure, and can be substan tial at microwave frequencies 
iLazarian fc Finkbeineill2003f) . 
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Finally, in the presence of a temperature quadrupole, 
the scattering of the CMB photons by free electrons in 
galaxy clusters leads to a secondary linear polarization. In 
addition to the primary CMB quadrupole, another impor- 
tant origin of temperature quadrupole is due to the pecu- 
liar velocity of clusters. However, the contribution to the 
B-mode polarization from this secondary signal is consid- 
erably low, several orders of magnitude below the lensing- 
induced polarizatio n on the large scales (for a review, see 
ICoorav et aT]l2004ll . 



3 GALACTIC FOREGROUNDS 

SUBTRACTION AND RESIDUAL POWER 
SPECTRUM 

Figure |5| clearly shows that the contribution of Galactic 
foregrounds to the B-mode polarization dominates over the 
CMB signal at all frequencies, even in the most optimistic 
cases. The knowledge of foreground emissions and the abil- 
ity to remove them assume therefore a primary role for the 
possibility of observing the polarization induced by grav- 
itational waves. Multifrequency observations are thus ab- 
solutely necessary, expecially on large scales. During last 
years, many methods have been suggested to perform the 
component separation of the microwave sky. Some of them 
adop t Bayesian approaches ( e.g., MEM by iHobaon gLaL 
I1993 and Wiener filtering by iTeemark fc Efstathioul Il996l 
and iBouchet fc Gis pert''l99u ) whereas others blind tech- 
niques ( e.g., SMICA by Dclabrouillo ct al. 2002aj and Fas- 
tlCA bv lMaino et"aD r2002.:'l. In this work we do not want 
to go into details of different techniques and assumptions 
for the components separation. Our aim is to estimate the 
level of residual Galactic foregrounds left in the clean maps 
after their subtraction. We make the following hypotheses: 
i) a template for synchrotron and dust emission (the only 
Galactic foregrounds considered) is available; ii) their spec- 
tral behaviour is estimated in the frequency range of in- 
terest. A linear subtraction of Galactic foregrounds is then 
employed extrapolating templates to the "cosmological" fre- 
quency. The instrumental noise present in templates and 
the uncertainty on the frequency spectrum of foregrounds 
are taken into account in order to compute the residual 
polarization. Compared to standard techniques for compo- 
nent separation, this method is, in principle, less accurate. 
However, it provides qualitative estimates for the residual 
foregrounds that are good enough for the aim of this work. 
Moreover, its simplicity allows to compute analytically the 
power spectrum of residual Galactic contaminants and, in 
addition, the method only requires minimal assumptions on 
the foregrounds and their statistical properties. 

Let Vo be the frequency chosen to detect the CMB polar- 
ization (hereafter called the "observational frequency") and 
ut the frequency of the template for the Galactic foreground 
we want to subtract. We suppose that the foreground fre- 
quency spectrum can be approximated by a power law (at 
least between Vt and Vo): 



1^0 



(5) 



In general, the value of the spectral index /3 is expected to 
change with the sky position. We indicate with 1^^ (n) and 



/3(n) the template foreground intensity and the measured 
spectral index in the sky direction n, while the difference 
between these quantities and the actual ones is A7^j (n) and 
A/3(n) respectively. If we extrapolate the template to the ob- 
servational frequency using Eq. the residual foreground 
left in the map at Vo after foreground subtraction is: 



A/,„(n) = J.„(n)-/,„(n) = 
(7.,+A/.,)(n)(^j 



-/3(n) 



Vt 



-A/3(n) 



- 1 



+ AJ.,(n)(^ 



-[/3{fl) + A/5(n)] 



(6) 



Keeping only first order terms, Eq. is reduced to 



AJ, 



.„(n)~ln(^)/.„(n)A/3(n) + A/.,(n)(^) . (7) 

As expected, the residual foreground in direction n consists 
in a part proportional to the uncertainty on the spectral in- 
dex, and another proportional to the error (or noise) on the 
foreground polarization at the template frequency. Extend- 
ing this result to the case of polarization measurements is 
trivial: 

A(0±i(/),„(n) ~ ln(^)(Q±if7),„(n)A/3(n) + 

-/3(n) 



A(Q±it/).,(n)p 



(8) 



Assuming that no correlation exists between the errors on /3 
and QiiziU and their estimated values, the X = E, B modes 
APS of the residual polarization is the sum of the two terms 
in Eq. The former (indicated by Cjj) is obtained by 
the convolution of the spectrum of the estimated foreground 
polarization at Vo, C{^^, with that of the uncertainty on the 
spectral index A/3, . The mathematical calculations are 
reported in Appendix A. We get 



C 



167r ^ 



^(2^i + l)Ci,^^(/,^i), 



(9) 



where A = Inii/t / Vo) and 



2 -2 / V2 -2 



(the last two terms are the Wigner 3j-Symbols). In prin- 
ciple, a convolution is required also for the second term of 
Eq. ||5J , C]^| . However, keeping again only first order terms 
(in particular, we replace /3(n) with the average value of the 
spectral index, {/3)) and supposing the polarization uncer- 
tainty in the template to behave like a white noise, we find 
that 



(5) 



-2(13) 



(10) 



where — 0,pixOt is the white-noise spectrum in the tem- 
plate map {^pix is the solid angle per pixel and at the noise 
per pixel in polarization). 

The hypothesis of a power-law frequency spectrum is 
not quite true for the dust emission. It can be better mod- 
elled by a law as Iv oc v'^ Bv{T), where Bv{T) is the Planck 
function at temperature T of dust grains. Then, for the dust 
emission we replace Eq. Q by 
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Figure 3. Residual polarization spectra under the hypoth- 

esis of average (CI; solid line) and pixel-dependent (C2; dashed 
line) spectral index for Galactic foregrounds. The dotted line is a 
power law normalized to the unity and with slope —1.8. 

We can assume the temperature of grains to be constant over 
all the sky and the spatial changes in the frequency spectrum 
to be restricted only to the value of /3. The residual power 
spectra for the dust emission will be still given by Eq. (j^Jl, 
but multiplied by the factor [(e'""/'^^ - l)/(e'"'°/-^'^ - 
that is ~ 2 for r = 18 K, Vo = 100 GHz and ut = 140 GHz. 



ror A/3(n) show coherent structures like those observed in 
total-intensity templates. This is visually observed compar- 
ing the maps o f the sy nchrotron spectral index provided by 
iGiardino et al.l ll2002l) and Berna rdi e t al. ( 2004) and the 
total-intensity ones. More quantitat ively, we estimate the 
power spectrum of /3(n) obtained bv IGiardino et all (|2002h 
from low-frequency surveys, and we find a power law spec- 
trum with a slope of ~ —3, in agreement with the behaviour 
of the synchrotron intensitv spectra ([T^ ^m^rk ^ Efstathioul 
llQQdiBouchet et alj[l993:lBouchet fc Gispertlll999l) . We ex- 
pect it is same for the dust emission. 

Therefore, in this case (CI) we assume that is pro- 
portional to the total— intensity APS for the Galactic 
foregrounds. The normalization of is obtained by fixing 
the dispersion of the spectral index around its average value, 
a/s, and using the relation: — ^^(2£-|-l)Cf Wf/47r, where 
We is the window function corresponding to the resolution 
at which the dispersion ap has been computed. 

The residual power spectrum obtained under this 
hypothesis is shown as solid line in Figure |3 Hereafter, we 
take: oc and aa = 0.2, in agre ement with the results 
found by IGiardino et^ (|2002h and iBernardi et"ai] t004\ 
for the synchrotron emission at high Galactic latitudes, and 
by WMAP me asuremen ts for the dust spectral index (see 
b igure 9 of iBennett et al.ll20oi^ . The polarization spectra 
of the Galactic foregrounds are assumed to be of the form 
^xe ~ Both spectra are smoothed by a window func- 

tion with FWHM= 1°. We notice that the residual spec- 
tra Cjl have the same slope of the foreground polarization 
spectra as expected (except for the first ^s) , but with an am- 
plitude a factor ~100 lower [here we have considered A = 1 
in Eq. ©]. 



3.1 Residual power spectrum due to spectral 
index uncertainties 

In order to estimate the power spectrum of the residual fore- 
ground polarization by Eq. linjl- llO^ , the knowledge of three 
angular spectra is required: (polarization spectrum of 
the Galactic foregrounds which can be directly computed 
from the template maps); (spectrum of the white noise 
in foreground templates); (spectrum of the spectral in- 
dex error). Here, we focus on the latter which is the most 
critical to define the behaviour of C^] . We consider two 
different extreme cases: 

(CI) Average spectral index: the template for a Galac- 
tic foreground can be extrapolated to the observational fre- 
quency using the sky-average spectral index. In this case, 
the error on the spectral index will be the difference be- 
tween /3(n) and the estimated average value {/3) = {(5) (as- 
suming no systematic effect in the estimation) . Because the 
spectral index of Galactic emissions is directly related to 
physical conditions of the region where the emission is pro- 
duced both for synchrotron (spectral index depends on the 
energy distribution of electrons) and for dust emission (de- 
pendent on the temperature and the emissivity of grains), 
we expect that maps of /9(n) and, consequently, of the er- 



(C2) Pixel dependent spectral index: as second case, 
we consider that the foreground template is extrapolated 
to the observational frequency using a spectral index esti- 
mated pixel by pixel. Moreover, we suppose that the error 
on these spectral indeces, A/3(fi), behaves like a white noise. 
Consequently its power spectrum will be constant. 
The amplitude of clearly depends on the noise in data 
used to compute the spectral index. As working case, we 
assume that the average error on spectral index is equal 
to the dispersion in the spectral index distribution, i.e. the 
above value cr/j = 0.2. Even if this assumption is probably 
very conservative, it is however interesting because it allows 
us a direct comparison to results of case CI, as the total 
power of the error is the same, but distributed differently 
in I. Figure |3 (dashed line) shows that the residual APS 
behaves like a white noise spectrum (at least up to the res- 
olution scale) with a very small amplitude at large angular 
scales: at ^ < 10 it is 10~^-10~^ times lower than C^^, and 
two orders of magnitude lower than predictions for case CI. 
Considering that the power of the CMB B-mode polariza- 
tion is concentrated at degree scales {I < 100), this second 
case for foreground subtraction is much more favourable, 
removing particularly well the Galactic foregrounds on the 
largest angular scales. 
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Figure 4. Polarization power spectra at 100 GHz (in tiiermodynamic temperature): CMB E— mode (upper solid line) and B-mode for 
r = 0.1 and 10~* (lower solid lines); synchrotron (left plot, dashed lines); dust (right plot, dashed lines); gravitational lensing (dotted 
lines); extragalactic radio sources with S < IJy (short dashed lines). The thick dashed lines show the residual contribution of the 
synchrotron and dust emission after their subtraction, following the two schemes described in sec. |3] In these spectra the effect of the 
beam (FWHM= 1°) start to be noticed at £ ~ 100. 



4 UNCERTAINTY ON THE 

TENSOR-TO-SCALAR PERTURBATIONS 
RATIO IN AN IDEAL EXPERIMENT 
INCLUDING FOREGROUNDS 

The main goal of CMB observations is to estimate the an- 
gular power spectrum. Its amplitude and shape allows us 
to have information on cosmological parameters. The pre- 
cision achieved for a parameter can be estimated from the 
Fisher matrix. The CMB polarization is usually described 
by E and B-mode and their relative spectra; for polarization 
measurements the cosmological parameter estimation using 
the Fisher matrix is given by 

( 1 dCEidCEi 1 dCBidCBi \ , . 
■^'^ \/\Cl, da, da, ^ AC|, da, daj ) ^ ' 

where ai are the cosmological parameters to be estimated. 
The diagonal terms of the inverse of the Fisher matrix pro- 
vide the minimum possible variance for the parameter ai, 
i.e. {SaiY ~ ^ii^ ■ The quantity ACxe is the uncertainty on 
the X — E, B power spectra at the multipole (.: 



{21 + l)fsky 



{Cxe + Nxef 



(13) 



where fsky is the fraction of the sky observed by the exper- 
iment. Beside the cosmic variance, the uncertainty on Cxe 
arises from the noise term Nxe, which includes the instru- 
mental noise and the foregrounds contamination, treated as 
an extra source of noise. 

As we have previously discussed, the CMB polarization 
at degree scales depends strongly on two cosmological pa- 
rameters, the tensor-to-scalar ratio r and the optical depth 
of the Universe r. Polarization spectra are affected also by 



other parameters, but polarization measurements are of par- 
ticular importance only for r and r, which can be at best 
poorly constrained with temperature. Therefore, we consider 
in the Fisher matrix only parameters at = {r, r}, all the 
other cosmological parameters being supposed to be already 
known. The results on the tensor-to-scalar ratio are given in 
terms of the quantity rHm, defined as the value of r under 
which the relative error on r is 0.3 (i.e., corresponding to 
about a 3-cr error) . On the contrary, the uncertainty on the 
parameter r will be not discussed, being its value generally 
estimated with a high precision, less than 10%, in almost all 
the cases we have considered. 

We consider as reference experiment a full-sky experi- 
ment with a resolution of FWHM=1° and a negligible level 
of instrumental noise ("noise-free"). A Galactic Plane cut 
is applied and fsky = 0.8 is used in the analysis. The obser- 
vational frequency, aimed to detect the gravitational wave 
polarization, is choosen at v =100 GHz, which corresponds 
to the frequency where the Galactic polarized emission is 
expected to be minimum (see Figure I^J. We model the po- 
larization of Galactic foregrounds using a power-law angu- 
lar spectrum: at 100 GHz, in antenna temperature, the syn- 
chrotron spectrum is 



Cxe = 1.2 X 10~^ r^-^ ifiK'^) 



(14) 



giving Xims = 0.12 /iK at 1° resolution (the normalization 
is based on the WMAP results and a polarization degree of 
20% ); for the dust polarization 



Cxe 



1.6 X 10~^^~ 



(15) 



cor responding to 0.26 uK a t 1° resolution (as expected by 
the iFinkbeiner et alJ 1^9^ model with a polarization de- 
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gree of 5%). These spectra, compared to the CMB E and 
B-mode ones, are plotted in Figure |1| 

We suppose that Galactic foregrounds are subtracted 
from data at 100 GHz using the prescriptions described in 
the previous section. In particular, we consider that a tem- 
plate for the synchrotron polarized emission is available at 
V = 60 GHz and one for the dust at 140 GHz. Under these 
conditions and supposing noise-free templates, the noise 
term in Eq. 11311 is given only by the residual Galactic po- 
larization, (thick dashed lines in Figure 4), and ex- 
tragalactic foregrounds. The rms amplitude of the Galactic 
polarization is reduced by a factor 20 or 60, according to the 
subtraction scheme used. 



4.1 Constraints on r imposed by Galactic 
foregrounds 

The Fisher matrix is used to estimate the uncertainty on 
r for our reference experiment. No contribution from extra- 
galactic sources is taken into account. The results are shown 
in Figure |S] The two different methods to subtract Galactic 
foregrounds are considered: in CI case (the frequency spec- 
trum for Galactic emissions is considered independently of 
the sky position) we find rum = 8.3 x 10~^; a detection of the 
tensor polarization with lower amplitude can be achieved if 
the spectral law is estimated pixel by pixel (C2 case). In 
particular, a 3-(t detection is possible for r = 1.8 x 10~^. 
These very low values for rum tell us that the Galactic fore- 
grounds are subtracted in a very efficient way in both cases 
(as expected when the noise is negligible). 

No improvement of the sensitivity on r is observed 
increasing the instrumental resolution to angular scales 
smaller than 1°, confirming that all the information on the 
cosmological B-mode comes from I ^ 100. If the resolution 
is reduced to FWHM= 7°, we find that rum gets clearly 
worse only in the CI case. On the contrary, the C2 results 
appear to be nearly independent of the resolution at least 
for FWHM< 7°, meaning that the largest scales {£ <, 30) 
are enough to provide a good detection of r. For a low- 
resolution experiment, the loss of information from £ ~ 100 
is probably balanced by a decrease of the amplitude of C^j . 

Finally, we have verified that previous results are nearly 
insensitive to the value of the optical depth. In fact, the 
reionization affects the B-mode APS only for £ ^ 10, where 
the cosmic variance strongly constrains the precision on the 
r detection. Hereafter, the results are given for r = 0.1. 

4.2 Effects of the extragalactic B— mode 
polarization on rum 

The next question is how much extragalactic radio sources 
and the polarization induced by gravitational weak lensing 
can affect the estimate of rum- At degree scales they provide 
only a secondary contribution to the B-mode polarization 
compared to the Galactic emissions. At 100 GHz, the radio- 
sources and lensing-induced B-mode spectra are constant 
with nearly the same amplitude, at least up to ^ ^ 300. How- 
ever, as shown in Figure|l] multifrequency experiments allow 
us to remove efficiently the Galactic foregrounds from po- 
larization maps, and thus extragalactic polarization becomes 
the greatest contribution to the B-mode also at very large 
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Figure 5. Fractional errors of r as function of r for the refer- 
ence experiment (full-sky, "noise-free" and FWHM= 1°) in the 
two hypothesis on C^: CI (solid line) and C2 (dot— dashed line). 
Dashed line (CI case) and dotted line (C2 case) are for the same 
experiment but with FWHM= 7° . Only Galactic foregrounds are 
considered. 

scales. Observations at different frequencies are not useful 
for the subtraction of point sources or of the lensing-induced 
polarization, that require experiments with arcminute res- 
olution. They become a critical problem for low-resolution 
experiments. 

4-2.1 Extragalactic radio sources 

In the last years, a large effort has been done to de- 
velop methods able to detect extragalactic point sources 
in C MB maps and to estimate their intensity emission 
(e.g.. iTeemark fc Oliveira-Costal Il998l . ISanz et al.l l200ll. 
IVielvaet al.ll2003^ . Complete catalogues of extragalactic 
sources with low flux limits are expected to be extracted 
from the forthcoming data. For example, IVielva et alJ 
l|20o3), using a method based on the Spherical Mexican 
Hat Wavelet, showed that 90%-completed catalogues up to 
~ 200 mjy can be extracted by Planck observations at the 
frequencies of 100, 143 and 217 GHz. At the moment, these 
techniques have not been extended yet to polarization maps. 

In any case, because we are only interested to reduce 
the contribution of point sources to the polarization maps, 
the best strategy is to exploit the total-intensity catalogues 
in order to mask all the pixels with bright sources. Suppos- 
ing that total-intensity catalogues of radio sources will be 
available up to very low flux limits, the only restriction is 
given by the number of pixels that can be masked in a map 
without losing too much information. This is a strong con- 
straint expecially for low-resolution experiments. We fix in 
the ten per cent the maximum percentage of beam-size pix- 
els we can mask. For a full-sky experiment with FWHM= 1° 
(like the reference experiment), it corresponds to 3300 pix- 
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Figure 6. Fractional error of r for the reference experiment, con- 
sidering the two different schemes of Galactic foregrounds sub- 
traction (CI, left plot; C2, right plot). The contribution of radio 
sources with S < 150 mjy (solid lines) and S < IJy is included 
(dashed lines). Previous results with no contribution from extra- 
galactic foregrounds are shown as the dotted line. 

els, taking into account the Galactic Piano cut. According 
to the number counts of Toffolatti ct al. (1998) model, the 
same number of radio sources is found at 100 GHz for fluxes 
5* > 150 mJy. 

Thus, we compute the relative error on r in the hy- 
pothesis that all the sources with S > 150 mJy have been 
masked (the APS amplitude of the remaining radio sources 
is 7.2 X 10"* nK^). We find ni„ = 1.4 x 10"" and 6 x 10"^ 
in the cases CI and C2 respectively (see Figure [SJ. In this 
figure we report also the results obtained if only the sources 
with 5 ^ IJy are removed. We stress that as the contribu- 
tion of extragalactic sources becomes higher the rum value 
is less affected by the amplitude of the residual Galactic 
polarization. 

Because the contribution of radio sources tends to de- 
crease at frequencies higher than 100 GHz, it is interesting 
to study what occurs when the observational frequency is 
increased, i^o > 100 GHz. We have considered i^o ~ 140 GHz 
and the Galactic foreground templates at 60 and 220 GHz. 
However, the resulting rum get significantly worse than in 
the reference case, confirming that 100 GHz is an appropri- 
ate frequency to detect the CMB B-mode polarization with 
this foreground subtraction scheme. 

4-2.2 Gravitational lensing 

The decontamination of B-mode polarization maps from 
the polarization induced by gravitational weak lensing can 
be achieved reconstructing the projected mass distribu- 
tion of la rge scale structures (or the lensing potential). 
IHu fc Ofa moto (2002) showed that the best estimator for 
the projected mass distribution is provided by the correla- 




r r 



Figure 7. Fractional error of r for the reference experiment, con- 
sidering the two different schemes of Galactic foregrounds sub- 
traction (CI, left plot; C2, right plot). We include the contribu- 
tion of the lensing— induced polarization: the total one (dashed 
lines) and reduced by a factor 10 (solid lines). No contribution 
from extragalactic foregrounds is considered as the dotted line. 



tion between the CMB E and B-modes, which allows us 
to reconstruct the distribution power spectrum up to mul- 
tipoles of I ~ 1000. Different a pproaches hav e been sug- 
gested by Knox & Sonej ^2002^ ■ iKesden et alJ ^2002^ and 
,Scliak & Hirata (2003): they show that a high precision in 
the lensing reconstruction can be obtained only for low-noise 
expe riments with arcminute resolution. In particular. Table 
1 of ISeliak &: Hiratal ll2003l) provides the residual B-mode 
contamination as a function of the beam size and of the in- 
strumental noise. The lensing "noise" always decreases when 
instrumental noise and/or the resolution are improved. The 
authors suggest that in an ideal condition of a free-noise 
experiment with unlimited resolution the lensing contami- 
nation could be removed completely. On the contrary, they 
argue that for wide beam detectors (FWHM> 20') the clean- 
ing methods are practically useless. 

So far, we have considered that the typical polariza- 
tion experiment is full-sky with degree resolution. In this 
case, the lensing contamination can become the dominant 
limitation for detecting the gravitational waves polarization 
with low values of r. A 3-o" detection, in fact, is possible 
only for r ^ 4 x 10"* if the lensing noise is not subtracted 
(see dashed lines in Figure |7|l. On the other hand, informa- 
tion from high-resolution experiments can be exploited: the 
Planck mission and upcoming ground-based experiments 
are expected to be able to measure the lensing-induced 
B-mode APS with a good precision (IBowden et al ] |2004l. 
[Canga 2004,). In the hypothesis that the lensing-induced 
spectrum can be red uced by a factor of 10 [as indicated by 
iKnox fc Sond (|20o3) for his reference experiment] , the pre- 
vious results are improved significantly, expecially in the C2 
hypothesis, for which we get rum — 5 x 10"^ (see Figure|7|l. 
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Figure 8. Fractional errors of r for the reference experiment, 
including all foregrounds: the two schemes for the Galactic polar- 
ization subtraction are considered (CI, solid line; C2, dot— dashed 
line). Dashed line (CI) and dotted line (C2) are for a small-area 
experiment without extragalactic contributions. 

4.3 Limits on r including all foregrounds 

Finally, we show in Figure|S|the results on Ar/r taking into 
account the following foreground contributions: 1) the resid- 
ual Galactic polarized emission; 2) the contribution from 
extragalactic radio sources with flux S < 150 mjy; 3) the 
polarization induced by weak lensing, but with its power 
spectrum reduced by a factor 10. For the reference exper- 
iment we obtain rum — 9 x 10~^ (1.8 x 10""*), using the 
C2 (CI) spectrum for the residual Galactic polarization. We 
want to stress again that, when extragalactic foregrounds 
are included, the value of rum is nearly independent of the 
method (and hypothesis) employed for removing the Galac- 
tic foregrounds polarization. If we were able to completely 
subtract the Galactic emission, the resulting rum (under 
previous conditions) would be 7 x 10~^, not far from the 
above values. It means that the limit of r ~ 10~* on the 
detectability of the polarization produced by gravitational 
waves is due essentially to radio sources and lensing, even 
in the case of a raw Galactic foregrounds subtraction (as in 
the CI scheme). 

When extragalactic contributions are included, low- 
resolution experiments are much less effective to estimate 
r. In particular, extragalactic radio sources are a serious 
problem due to the difficulty to detect and subtract them. 
For example, considering FWHM=7°, we find rum ~ 10""^ 
in the best case. 



^.3.1 High-resolution, small-area experiment 

An experiment with arcminute resolution could overcome 
the constraints imposed by extragalactic foregrounds. In 
fact, on the one hand, the number of arcminute pixels that 



can be masked to remove radio sources is much higher, re- 
ducing the flux limit to <C 150 mJy and making their contri- 
bution negligible (also in a small area); on the other hand, 
the high resolution (with the "noise-free" condition) allows 
us to reconstruct with high precision the lensing B-mode 
power spectrum directly from the data llSeliak fc Hiratal 
|2p03). In particular, we study what occurs for an experi- 
ment covering a sky area of 30° x 30° in the hypothesis of a 
perfect subtraction of extragalactic foregrounds: the limited 
sky area strongly constrains the detectability of the tensor 
polarization, due to the cosmic variance on the low I (we 
remind that the uncertainty on r increases proportionally 
to fJi^y'^)- Only if Galactic foregrounds are removed in an 
accurate way, as under the C2 conditions, rum is similar to 
that found for a full-sky experiment (see the dotted line in 
Figure |HJ . 



5 RESULTS FROM FUTURE EXPERIMENTS 

In this section we apply the previous analysis to experi- 
ments characterized by a given sensitivity to the polarized 
emission. Practically, it means to include the instrumental 
sensitivity in the noise term of the Fisher matrix. We also 
take into account the noise in template maps in order to 
estimate the APS of the residual Galactic polarization, by 
adding the term given in Eq. llOt . 

In the next years, a large number of experiments has 
been planned with the explicit objective to measure the 
CMB polarization (for a review on the polarization exper- 
iments sec Dclabrouillc ct al. 2003b). However, while they 
may provide an accurate estimate of the E-mode power 
spectrum both at large and small angular scales, for the 
forthcoming experiments the detection of the gravitational 
waves polarization is still a challenge and it will be possible 
only for high values of the tensor-to-scalar ratio r. 

In the mid term, the Planck mission* satisfies the best 
conditions to search the "inflationary" polarization: a full- 
sky coverage and a good characterization of all the relevant 
foregrounds. In fact, a large frequency range is spanned and 
more than one frequency is dedicated to the study of the 
Galactic foregrounds, providing adequate knowledge of their 
spectral behaviour and template maps in polarization. On 
the other hand, the arcminutes resolution allows to detect 
point sources up to quite low fluxes and to have a first esti- 
mate of the gravitational lensing contribution. Nevertheless, 
the instrumental noise is still a very strong constraint: the 
sensitivity to Stokes parameters Q and U per square de- 
gree pixel is higher than 1 /iK at all the frequencies (after 14 
months) , significantly higher than the expected rms value for 
the cosmological B-mode polarization. Taking the 100 GHz 
channel as the "observational frequency" (and using the 30 
and 217 GHz channels as templates for the synchrotron and 
dust emission respectivelly) , we find that a 2-a detection of 
the gravitational waves polarization is possible for r ~ 0.05. 
This result does not change if we do not include the contri- 
bution of foregrounds. 

The research of the gravitational waves polarization 



* http:/ /www. astro. esa.int/SA-general/Projects/Planck/ 



12 Tucci et al. 




(7/a{Planck) 

Figure 9. The value of r;i,„ for an experiment lilce Planck, but 
with a sensitivity improved by a factor o-/apia„ck- In the upper 
(lower) line Galactic foregrounds are subtracted following the CI 
(C2) scheme and extragalactic foregrounds are partially (com- 
pletely) removed. 

is also promising using ground-based experiments. As dis- 
cussed above, the disadvantage of covering a limited area of 
the sky can be partially compensated by the high resolution, 
which permits a better control of the contribution provided 
by extragalactic sources and gravitat ional len sing. Experi- 
ment s like PolajB ear (see | Gangal20o3) , Clover jTavlor et alJ 
|2004) or QUIET jMilleill2004) are planned to achieve a very 
good sensitivity, of the order of 0.1 /iK or less, using arrays 
of 1000s of detectors in only two or three frequency channels. 
As example, we discuss in detail the case of PolarBear: it will 
operate at the frequencies 90, 150 and 220 GHz covering 225 
square degrees; the number of detectors for each frequency 
should be 400, 600 and 200 with a sensitivity per detector 
of 310, 345 and 640 ^Ks^''^ respectively. Considering an ob- 
servational time of 1.5 years with an efficiency of 30%, the 
sensitivity per 1° pixel at 90 GHz would be ~ 0.1 ^K, an 
order of magnitude lower than the Planck one. We find that 
PolarBear will be able to provide a 2-a detection of the 
tensor polarization for r ~ 0.07. 

Based on the above experiments, we describe below the 
possible features for future missions specifically dedicated 
to the detection of the CMB B-mode polarization, and we 
compute the corresponding detectability limit on r. 

Space experiment: it is interesting to know how the de- 
tectability limit on r improves when reducing the instru- 
mental noise in the Planck experiment by the same factor 
at all the frequencies. Figure |3| shows the results for rum - In 
the upper line we are supposing that Galactic foregrounds 
are subtracted following the CI scheme and extragalactic 
foregrounds are partially reduced (as discussed in Sec. 4.3). 
This is probably an upper limit. On the other hand, in the 
lower line extragalactic foregrounds are neglected (we con- 
sider they are completely removed), while the C2 scheme is 



used for the Galactic polarization subtraction. We observe 
that the results are independent of the two hypothesis on 
the foregrounds subtraction if Planck-detectors sensitivity is 
improved by a factor ^ 30 (corresponding to rum ^ 10"'^). 
In this case, the instrumental noise is still the most con- 
straining contribution. For sensitivities 100 times or more 
better than the Planck one, the upper line is constant with 
riim — 3.5 X 10~*, while the lower line still decreases reaching 
a riim similar to the free-noise case for a/apianck ~ 10"''. 
We plot the tail of the lower curve as a dotted line in or- 
der to indicate that, for so high sensitivities, it must not be 
considered a lower limit for the detectability of r, since the 
uncertainty of foregrounds spectral index we used, = 0.2, 
can be improved. 

We discuss now a more concrete example for a space 
experiment: a total of 10,000 detectors in three channels at 
100, 140 and 220 GHz with the same sensitivity of Planck 
detectors. With about 3,000 detectors for each channel, in- 
stead of 8 present in Planck, we expect the sensitivity to 
increase by a factor of about 20. The 100 GHz channel is 
used as observational frequency and the other two channels 
to remove the dust polarized contribution. A template for 
the synchrotron polarization is also required, being its signal 
still important at 100 GHz. If the low-frequency channels of 
the Planck mission are considered, we find rurn — 2.5 x 10""'. 
In spite of the low accuracy in the synchrotron subtraction, 
riim is only a factor ~ 1.4 higher than the corresponding 
value of Figure |U] The results are not affected by the contri- 
bution of extragalactic foregrounds (radio sources and lens- 
ing) , while only a slight dependence is noticed on the scheme 
of the Galactic foregrounds subtraction. 

Ground Based Bolometers: an array of 5,000 bolometers 
at the frequencies 90 and 150 GHz with a sensitivity per de- 
tector of 250 /iKs^''^. After one year of observations (60% 
efficiency) over a 30° x 30° patch of the sky, we obtain a 
sensitivity of 33 nK per 1-degree pixel at both frequencies. 
Using the 30-GHz channel^ of Planck mission to subtract 
the synchrotron emission, we find rnm ~ 0.03. These val- 
ues can be reduced by a factor 2 if a channel at 100 GHz is 
choosen instead of 90 GHz. The results are practically inde- 
pendent of the scheme to subtract Galactic foregrounds and 
of the contribution of extragalactic foregrounds. Reducing 
the area observed, the sensitivity per pixel is improved but 
not the sensitivity on r: for a 20° x 20° area, a = 22 nK 
while riim — 0.07. 

Ground Based Radiometers: an array of 1000 radiome- 
ters at 30 GHz and 4000 at 90 GHz. Assuming 1 year of ob- 
servations over a 30° x 30° patch with a 80% efficiency and a 
sensitivity of 120 and 200 fiKs^^^ at 30 and 90 GHz respectiv- 
elly, we get sensitivities per square degree pixel of 31.5 nK at 
30 GHz and 21 nK at 90 GHz. Using the 217-GHz channel'^ 
of Planck for the subtraction of the dust emission, we find 
a 2-a detection for r ~ 0.28. This result points out that 
experiments based only on radiometers, although with com- 
parable or better sensitivity, appear to be less efficient than 



^ In this case, the residual polarization arises mostly from the 
noise in the template map, rather than from the uncertainty on 
the spectral index. So, foreground templates at frequencies distant 
from the observational one are preferred. 
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bolometers to detect the tensor polarization [or, in other 
terms, that higher frequencies {v ^ 100 GHz) are better 
than lower ones (z/ ^ 100 GHz)]. The reason is the different 
spectral behaviour between dust and synchrotron emission: 
the frequency spectrum of synchrotron is very steep and al- 
ready at 100 GHz its polarized contribution is small; on the 
contrary, the dust spectrum decreases quite slowly with the 
wavelength and at 70-80 GHz its polarization level is still 
comparable to the synchrotron one. 

Finally, considering jointly radiometers that operate at 
30 and 70 GHz and bolometers at 100 and 150 GHz with the 
features previously described, we find a sensitivity to the 
CMB B-mode polarization similar to the one achieved by 
the space mission: rum ~ 6 x 10~^. 



6 DISCUSSIONS AND CONCLUSIONS 

In this paper we have investigated the constraints that fore- 
grounds can impose on the detection of the B-mode po- 
larization induced by gravitational waves. First of all, we 
carry out a detailed analysis on the expected level of the 
foregrounds polarization at degree scales, which are the an- 
gular scales where the cosmological B-mode signal is con- 
centrated. The diffuse Galactic emissions (synchrotron and 
dust emission) are the dominant components. The mini- 
mum of their polarized emission is found to be around 
100 GHz, corresponding to a rms value for the B-mode of 
~ 0.5 jUK at 1° resolution (to be compared with the CMB 
i?rma — 0.43-^ /iK) . We expect a shift of the "cosmological 
window" in polarization towards higher frequencies as com- 
pared to temperature fluctuations, where the lowest emissiv- 
ity of foregrounds is observed around 60 GHz. It is due to the 
high polarization degree of the synchrotron emission, whose 
diffuse component is observed to be up to 30% polarized or 
more. The contribution of extragalactic radio sources and 
of the polarization produced by gravitational lensing is sec- 
ondary on large scales, but not negligible when techniques 
to subtract Galactic foregrounds are applied to polarization 
data. On the contrary, other sources of polarized emission 
(such as infrared galaxies, SZ effect, etc.) are expected to be 
negligible. 

The subtraction of the Galactic contributions from po- 
larization maps is mandatory to search for the gravitational 
waves signal. Assuming component separation boils down to 
subtracting, at the CMB observing frequency, a set of fore- 
ground templates scaled in frequency according to a model 
emission law, we provide an analytical expression for the 
power spectrum of the residual polarization, depending on 
the noise in the foreground templates and on the uncertainty 
in their spectral behaviour. We study two concrete cases: i- 
the frequency spectrum of Galactic emissions is supposed 
independent of the sky direction; ii- the frequency spec- 
trum is estimated pixel by pixel, but with an uncertainty 
that behaves like white noise. As expected, under this sec- 
ond hypothesis a better subtraction of Galactic foregrounds 
is achieved: if only the contribution of the residual Galactic 
foregrounds is taken into account, we find that a 3-a de- 
tection of the tensor polarization is possible for r ~ 10~^, 
corresponding to an energy of inflation of ~ 2 x 10^^ GeV. 
However, this value should not be considered as a lower limit 
for the detectability of the primordial B-mode polarization 



since it arises from working hypothesis. In fact, in this case, 
the amplitude of the residual Galactic APS depends on the 
average uncertainty with which the spectral index is com- 
puted. We have taken cr^ — 0.2: this value is probably reli- 
able for the data available up to now, but it will be improved 
in the future. In this sense, a limit on the detectability of 
the CMB B-mode polarization cannot be put "a priori" . Op- 
posite conclusions are found if the subtraction of Galactic 
foregrounds is based on the assumption that their spectral 
law is independent of the sky direction. Now, the amplitude 
of the residual Galactic foregrounds depends on the intrinsic 
dispersion of the spectral law in the sky (cr^ — 0.2) and not 
on our ability to estimate it. So, the value of r ~ lO"** that 
we find for a 3-cr detection corresponds to the lower limit 
for the detectability of the gravitational waves polarizationin 
this case. 

Clearly then, taking into account the spectral variabil- 
ity of foregrounds in the component separation is a key 
point in order to detect cosmological CMB B-mode polariza- 
tion for very low r. At the present time, only few methods 
have been proposed to perform the component separation 
with spatial variabil ity of foregrounds dStolvarov et al .120041 : 
iBennett et"ai]|2003l'l showing promising results. More efforts 
within this field are still required. 

The effects of extragalactic radio sources and gravita- 
tional lensing on the r sensitivity is also studied. After sub- 
traction of the Galactic foregrounds, they become the major 
contaminant also on the large scales. In principle, CMB ex- 
periments with high resolution can control accurately these 
foregrounds. Lensing and radio sources are much more prob- 
lematic for measurements with a wide beam and comple- 
mentary data with high resolution are required to reduce 
their polarized contribution. We have seen that, for a full- 
sky "noise-free" experiment with 1° resolution, a reduction 
by a factor 10 for the B-mode APS of radio sources and 
lensing is probably feasible, giving a detectability limit of 
r ~ 10"'' independent of the scheme to subtract Galactic 
foregrounds. This value is, in any case, better than those 
found by experiments covering only a small area of the sky, 
that are strongly limited by the cosmic variance. 

Finally, we have discussed what are the prospects, "in 
practice" , for the detection of tensor polarization in the fu- 
ture. In the near term, a tensor-to-scalar ratio of r ~ 0.05 is 
well within reach of the Planck mission and planned ground- 
based experiments. A detection of tensor modes with so high 
values of r would be an indication of an infiation driven by 
some exotic physics llKinnevll2003ft . The next generation of 
missions, that are at the moment under discussion by NASA 
and ESA, will be able to significantly improve the sensitivity 
on the polarized emission by means of arrays of 10^-10'' de- 
tectors, and will probably allow us to investigate more con- 
ventional theories of the infiation, spanning ranges of energy 
down to ~ 5 X lO'^GeV. 
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APPENDIX A: ESTIMATE OF THE RESIDUAL POWER SPECTRUM 

In this appendix we demonstrate the result given in Eq. @. In particular, we want to calculate the E- and B-mode power 
spectrum for the function 7?.{n) = 7'(n)/3(n), where P is a spin-2 function [in Eq. Q it is the linear combination of Stokes 
parameters Q ± iU] and /3 is a real scalar function (i.e., the error on the spectral index). 

A spin-2 function can be expanded on t he sphere by the spin-2 sp herical harmonics: /(n) = X/fm '^^/m 2Yirn,(.'n) and 
fin) = X/£m '^-^.fm -2y«m(n). Accordiug to lZaldarriaga fc SeliakI (|l993) formalism, the _B-mode power spectrum is 

2 1 , * 

CbI =< \aBlm,\ > ~ ^ < {a,-2,£m — 0.2,lm)(a -2,lm ~ 0,2,lm) > = 

= 4( 0,2,lma2,lm > + < ^-2/ m<ll2,f m > — < a-2,«ma2/m > — < 0.2,lma*-2,lm > ^ ■ (^1) 

In order to compute C^i, we start by considering the first term in the right hand of Eq. lAlll 

<ali.maf}m >=< / dr!2y/™(n)P(n)/3(n) / d^i' _2y/-™(fi')^*(n )/3(n ) > ■ (A2) 



Using the appropiate expansions for ^(n) — X/£m "^Mm 2^im(n) and /3(n) = X^fm ^fm^*'"(^)' ^'^'^ removing from integrals 
all the terms independent of the sky position, we find 



< af^emaf'im > = < X/ X/ "J^imi^fama J dQ2Yi^{n)2Ye^jni{n)Ye2m2in-) ^ ^ 

£imi £21^2 i^m^ £47714 

X J dn' -2Y/_^{n') -2Ye^,n.^{n)Yi^miin) >= X! X! X! X! < 



v p 



-2,f3m3 >< a-i^^^a^^^^ > X 



l\m-i ^2"^2 ■^3''^3 ^4^JT'4 

rfr2 2y«m(n) 2y(?imi(n)y«2m2(n) / df^' -2y«*_m(n') -25^f3m3(n')Y«^„4(n') . (A3) 



Now, taking into account that 2Yi^ = _2y/_„ , < a^fimi"-2,£3m3 >= (^)'"^ < laJ^imiP > <5£if35,„i_m3 and 

< "Im2«lm4 ^f/^ -rn^ , We obtain 

, Tl TZ* , _ \ ^ \ V |2 ^ r^l3 ^ 

< "2,fmffl2,^m > — / ^ < P2/imil > ^ 

X J^i-r'^""^ J dn.2Ye^^{n)2Ye,^,{n)Y,^,in) J dfi' 2y,™(n') -2y«i-mi (fi')V^2-™2 (" ) ■ (A4) 



mi m2 

According to properties of spin-weighted spherical harmonics, we have 



V47r L J V— si — S2 — S3/ VTii m2 ms/ 



The last two terms are the Wigner 3j-Symbols and verify the condition 

l2 

Vmi 7712 ma/ 



■^1 ^2 ^3 \ , „ .J. J mi + m2 + ma = 
if 



and the orthogonality relation 

V(2.. + l)(^^ '^)('^ ^M=l. (A6) 
^ — ' \mi m2 m-iJ \m\ m2 ma/ 

"127713 

Hence, Eq. (IA4II can be simplified 

< '^2,^m^i2/m > — / ^ < !^i2,£imij > ^ / ^ \ — ) X 

El£2 mim2 

^ V2 -2 / V-m mi m2/ V-2 2 /Vm-mi-m2/ 

sEp'. + ■) < I"" f > 12'"' + % ) ("2 ) ■ 

£1 I' 



It is easy to verify that an equivalent result is found for all the other terms in Eq. ijAlfl . Finally, we get: 

We remind that the sums on £i and are done only for values j£i — ^2] ^ i' ^ i'l + ^2. The same result can be found for the 
E-mode spectrum. 



